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Instantons are the natural mechanism in non-perturbative QCD to remove helicity from valence quarks and 
transfer it to gluons and quark-antiquark pairs. To understand the extent to which instantons explain the so-called 
"spin crisis" in the nucleon, we calculate moments of spin-dependent structure functions in quenched QCD and 
compare them with the results obtained with cooled configurations from which essentially all gluon contributions 
except instantons have been removed. Preliminary results are presented. 
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1. MOTIVATION 

In recent years there has been a major ex- 
perimental effort to measure structure functions 
which characterize the distribution of quarks in 
the nucleon. The challenge for theorists is now 
to understand the data which already exist and 
to predict the structure functions which will be 
measured in currently planned experiments. Lat- 
tice QCD provides the only known framework 
for non-perturbative calculation of hadron struc- 
ture, and although it cannot address the intrin- 
sically Minkowski structure functions themselves, 
through the operator product expansion it is pos- 
sible to calculate their moments. There have been 
several calculations of these moments in recent 
years, including calculation of all moments of the 
spin independent and longitudinal spin depen- 
dent structure functions through order four [Q 
and the tensor charge ^ by the QCDSF collabo- 
ration, as well as calculation of connected and dis- 
connected contributions to the axial charge H|J 
and to the tensor charge ||. These results agree 
qualitatively with experiment and discrepancies 
may plausibly be attributed to a combination of 
finite lattice size effects, and the omission of sea 
quarks and disconnected diagrams. 

However, in addition to showing that numeri- 
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cal solution of QCD reproduces the experimental 
results, we want to understand as fully as pos- 
sible the physical origin of the observed struc- 
ture. Since instantons play a major role in the 
physics of light quarks dJ?]], we also would like 
to understand their role in hadron structure func- 
tions. Especially interesting from this perspective 
is the spin structure of a proton. The 't Hooft 
instanton interaction UBULdndLSRSL is the only 
known vertex in QCD that directly removes helic- 
ity from valence quarks and transfers it to gluons 
and quark-antiquark pairs, thereby rendering in- 
stantons the natural mechanism to explain the 
so-called "spin crisis" . 

Hence it is important to identify the instanton 
contribution to the proton spin structure func- 
tions. On the lattice, a practical way of extracting 
the instanton content is "cooling" [|||] , in which 
one sequentially minimizes the action locally on 
each link and iteratively approaches a stationary 
solution. 

In this work we calculate the zeroth moments of 
the <7i(x) and h\(x) structure functions with and 
without cooling for u and d quarks in quenched 
lattice QCD. The zeroth moment of gi(x) is pro- 
portional to the quark's spin contribution to the 
total spin of a nucleon in the parton model and 
has been measured experimentally. The h±(x) 
structure function characterizes the transversity 
of quarks in a protonM. 
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2. METHOD 

We use a lattice of size 16 3 x 32 with f3 = 6.0 
and generate quenched ST/ (3) gauge fields using 
overrelaxed heatbath Monte-Carlo with Cabibbo- 
Marinari decomposition. We calculate the follow- 
ing matrix elements: 

(PS\qY*J 5 q\PS) = 2S i Aq (1) 

(PS\qa 0i ij 5 q\PS) = 2S l 6q (2) 

where IPS') is a ground state of a proton with 
momentum P and spin S (P 2 = Mp and S 2 = 
—Mp), and q defines the quark flavor u or d. 

The axial charge is related to g\ (x) through the 
operator product expansion by 

2 [ dxg 1 (x,Q 2 ) = 
Jo 

]T e[% 2 /Q 2 , g(p)) Z A ((atf,g(a)) Aq(a) 

q—u : d 

where e\jo is the Wilson coefficient and Za is 
the lattice renormalization constant for the axial 
charge and analogously for the tensor charge. To 
compare below with experimental data and with 
the calculation of rcf.|l|, we use the perturbative 
value of = 1) = 0.867 used in their work. 

The nucleon matrix elements are extracted 
from the ratio 

m = r-<«,> (3) 

where 0(x, t) is the operator of interest, either 
0{x, t) — qYi"j5q(x,t) or the tensor operator; the 
spin vector is chosen to be S = (0,0,0, |S|); the 
corresponding spin projection operator is T = 
(1 + 7 5 7 3 )/2; the nucleon source is J a (x) = 
£ab c Ua(u b Cj5U c ); and the states are normalized 
by (P\P) = 2EV6p p, 

For the hadron source (J) we use gaussian 
smearing with Vr 2 » 0.5 fm in combination with 
Coulomb gauge fixing at the time slice of the 
source. We use a point-like sink (J') and because 
there is no need for finite P to avoid operator 
mixing, we project onto zero spatial momentum 
to minimize numerical noise. 

We evaluate the numerator in (^|) using a stan- 
dard sequential source to create a set of back- 
ward propagators which can be contracted with a 



given operator and smeared forward propagators 
to produce the corresponding matrix element. 

3. CALCULATION 

For uncooled QCD we use k = 0.1515 and 
k = 0.1550, where for reference we note that for 
(3 = 6.0 k c = 0.1569 We define the lat- 

tice spacing from the SESAM nucleon mass chiral 
extrapolation yielding a^ 1 — 1.95 GeV. At this 
value of a our k's correspond to w 590 MeV 
and m-H w 980 MeV respectively. 

In the cooled case we use 25 cooling steps as in 
ref. Standard chiral extrapolation yields the 
following results: 



Table 3 



Chiral extrapolation in the cooled case 



K 


am n 


arriN 


0.1230 


0.405(12) 


0.748(24) 


0.1250 


0.278(17) 


0.615(31) 


k c =0.1269 




0.497(38) 



Using arriN to fix the lattice spacing we obtain 
a" 1 = 2.2 GeV. 

The three point function was evaluated with 
the source at t = 10 and the sink at t = 24, 
and the operator was averaged over the central 
plateau from t = 12 to t = 21. The statistical 
errors in R(t) in Eq.(||) were determined by the 
jackknifc method. Because the initial sample size 
for the uncooled configuration was only 40, lead- 
ing to some uncertainty in the optimum size of 
the central plateau, to be conservative we added 
a 5% systematic error to the purely statistical er- 
ror. 

4. RESULTS 

Our results for the unrenormalized axial vec- 
tor charge (An, Ad) and tensor charge (Su, Sd) 
are shown in Table 1 for uncooled QCD and in 
Table 2 for cooled QCD. 

In order to compare with QCDSF calculations 
and with phcnomcnological data, we consider the 
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Table 1 



Unrenormalized axial vector charge (Am, Ad) and tensor charge (Su, Sd) for the un-cooled case 



K 


Am 


Ad 


Su 


Sd 


0.1515 


1.118(43) 


-0.281(19) 


1.297(40) 


-0.302(26) 


0.1550 


0.967(124) 


-0.318(70) 


1.210(225) 


-0.329(58) 


Table 2 

Unrenormalized axial vector charge (Am, Ad) and tensor charge (Su, Sd) for the cooled case 


K 


Am 


Ad 


Su 


Sd 


0.1230 


0.814(22) 


-0.230(10) 


0.969(19) 


-0.236(10) 


0.1250 


0.703(48) 


-0.232(27) 


0.916(42) 


-0.222(30) 
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Figure 1. Lattice results for the axial vector 
charges 

axial vector charge, for which the renormaliza- 
tion factor Zji(l,g = 1) = 0.867 is available, 
and show our results in Figure 1. For uncooled 
QCD there is excellent agreement between our 
renormalized results for ZaAu and Z^Ad and 
the high statistics QCDSF calculation. Both 
lattice calculations extrapolate well to the phe- 
nomcnological data, indicating that at least for 
these moments, quenched QCD is qualitatively 
adequate. To understand how much of the axial 
charge arises from instantons alone, we also show 
the cooled results Am and Ac? in the same fig- 
ure. In principle, one should calculate renormal- 
ization factors arising from instantons. However, 
although we have not yet done this, we believe 
that most of the renormalization arises from fluc- 
tuations which have been removed by cooling, so 
that the relevant renormalization factor is close 
to unity. Hence, we believe that the qualitative 



agreement between the cooled and uncooled re- 
sult in Figure 1 strongly suggests that instantons 
indeed provide the dominant contribution to the 
nucleon spin. 
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